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ABSTRACT 

A model of the soft X-ray background (0.1 & hv < 0.284 keV). is 
presented in which the Sun is assumed to be inside an active supernova 
blast wave. The blast wave, in turn, is evolving in a preexisting 
cavity. The broad-band surface brightnesses can be explained by such 
a blast wave with an explosion energy of E Q 5 x 10 5 ® ergs and radius 
80 to 100 pc, using 6olar abundances. 

An approach to treating the problem of large anisotropies in the 
ambient medium is also explored, accommodating the observed 
anticorrelation between the soft X-ray surface brightness and the 
21 cm column density. It is found that only for post shock 
temperatures below 10^ K a shock propagating into a density 
enhancement will be dimmer than a similar shock in a lower density 
region. 

I. Introduction 

* 

This work continues a project of modeling the soft X-ray 
background by the assumption that the Sun is inside an active blast 
wave. Cox and Anderson (1982; hereafter CA) produced such models in a 
uniform ambient medium, and found a set of dynamical parameters which 
gives a tolerable fit to the sky surveys of McCammon et al. (1983) and 
others. The present work makes an attempt to generalize the CA 
project, finding limits on the parameters in question. Using 
dynamical models of Cox and Edgar (1983; Paper I) and Edgar and Cox 



(1984; Paper II), it is possible to examine the effects of a 
pre-existing cavity (say, from a previous supernova in the same 
vicinity) on a supernova blast wave. 

The McCammon et. al. (1983) sky survey is presented in 
broad-band form. The two lowest energy bands, Boron (B-band; 
0.1 < hv < 0.187 keV) and Carbon (C-band; 0.1 < hv <0.284 keV) , can be 
fit well by a model such as that discussed above. However, CA were 
unable to create a model which produces sufficient medium-energy 
(M-band; 0.4 < hv < 1.4 keV) counts. While some substantial fraction 
of the observed M-band sky brightness may come from sources other than 
hot Interstellar gas (for example dll stars or some extragalactic 
source; McCammon et. al. 1983), the M-band sky is remarkably 

isotropic outside of the Loop I/North Polar Spur complex, which 
suggests some local isotropic source such as a blast wave. Also, 
since the North Polar Spur is proposed (Borken and Iwan 1977, Iwan 
1980) to be a re-heated supernova remnant (SNR), and it is clearly 
associated with an M-band enhancement, it is desirable to create a 
model which will produce more M-band flux than those of CA. Toward 
this end, and also with an eye toward examining reheated cavity 
models, the present work has focused mainly on blast waves propagating 
into pre-existing cavities. However, using the solar abundances of 
Ross and Aller (1976), our models are no more productive of M-band 
photons than those of CA. 
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Section XI discusses the detailed numerical models and presents 
an analytical approximation to the predicted surface brightnesses in 
the B and C bands. Section III is an exposition of a scheme for 
examining a non-isotropic ambient medium in an attempt to model the 
observed anticorrelation between soft X-ray intensity and neutral 
hydogen column density. In section IV the results are summarized in a 
series of contour plots, and limits on the dynamical parameters are 


discussed 
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II. Detailed Numerical Modeling 

We have generated a modeling program which follows the dynamics 
of a blast wave using the dynamical model of Paper II, integrates the 
ionisation balance equations using rates given by Raymond and Smith 
(1977, 1979), and then calculates the X-ray surface brightness. This 
is calculated in the form of binned spectra and integrated against 
band response functions, in particular those appropriate to the 
Wisconsin rocket detectors, as given by McCammon et al. (1983). We 
report here the band-integrated spectra, which facilitate comparison 
to existing sky survey data. 

We have also Included the preliminary response function for a 
Wisconsin rocket payload using beryllium filters to delimit the 
softest energy band (0.07 < hv < 0.110 keV). Preliminary results and 
a brief description of this payload are given by Sanders et 
al. (1984). 

This model assumes that the post-shock ion and electron 
temperatures are equal (electrons being heated by non-Coulomb 
processes In the shock). Early In the life of the remnant, thermal 
conduction is important, and is accounted for as recommended by Cowie 
(1977), allowing for saturation effects. As the post-shock 
temperature drops with time, the thermal conductivity also drops, so 
the structure of the model remnant approaches that of an adiabatic 
blast wave (see, e.g. Sedov 1959, and Cox and Franco 1981). 
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The characteristic radius for the changeover from early 
conduction dominated behaviour (described by Paper I) to late 
adiabatic behaviour is called R^. See Paper II for its prescription. 
It is frequently useful to parameterize the shock radius in units of 
Rj^ introducing z«R 8 /R 1 (E 0 ^ ,nj ) . (The algebra is complicated in the 
case of a pre-existing cavity by the need to specify the density at 

f 

some particular radius. It is most convenient to specify n^, the 
ambient density at R B - R^, putting one in the odd position of 
calculating R^ in terms of one of the parameters there. 

The models assume a power-law ambient density, considering 
n. 0 o R“«, with to ■ 0, -2, and -4 (uniform density and two steepnesses 
of pre-existing cavities). Solar abundances were used, and pre-shock 
ionization states were those of an equilibrium 1.0 x 10^ K gas. 

We find several interesting results, which tend to confirm those 
of CA. Perhaps most notable is that the various ions studied (all 
stages of the elements He, B, C, N, 0, Ne, Mg, SI, S, Ca, Fe and Ni) 
are typically far from colllsional equilibrium; time scales for 
ionization and recombination are comparable to the dynamical time 
scales. We also note that the ambient density is the most tightly 
constrained parameter, though the particular value one must assume to 
produce a desired surface brightness depends somewhat on the value 
assumed for u). 

Table 1 gives the dynamical properties of the runs that were 
done; the X-ray surface brightnesses (as viewed from the center) are 
given in Table 2. Nearly all the’ runs were for w ■ -4 since the 
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primary purpose of this effort was to examine the deviations from the 
CA calculations. 

We have found an analytic approximation to both the C-band rates 
and the B/C ratio, which is useful for discussing the behavior of the 
results. The C-band approximation is a slight modification of the 
expressions given by Hamilton, Sarazin, and Chevalier (1983; HSC) and 
CA for the total luminosity of a blast wave (with x s - as the 
density compression factor at the shock): 
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with a further factor of exp (1-10^K/T S ) when T s < 10^ K. This last 
factor (the "cold correction") is an attempt t!o reproduce the fact 
that lower temperature plasmas tend to radiate in lower energy 
portions of the spectrum than the C-band. The various factors of this 
formula can be given physical interpretations as follows. The 
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(where Tg - T Q /10 6 K) law 18 the equilibrium cooling coefficient of 
Raymond, Cox, and Smith (1976), appropriate when ionization stages are 
in approximate collisional equilibrium. The other alternative is the 
estimate by CA of a non-equilibrium cooling coefficient, modified by 
an exponential factor based on that of HSC which corrects for the 
finite time needed to reach ionization stages which will radiate soft 
X-rays. The oj-dependence was introduced to fit the w ■ -4 results of 
this paper, while preserving the to *■ 0 form of HSC. Further, an 
empirical fit to the B/C ratio produces 


5. « min [0.58, (- 
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-0.4386 r gQ cm “3 y r ^ 
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( 2 . 2 ) 


Since most of the runs in this work were done with to ” -4, the 
(o-dependence of equation (2.2) should be taken as uncertain. 

The results of these approximations as applied to the detailed 
runs are presented in Table 2. Also included are the percentage 
differences between the numerical models and the analytical 
approximations. While the agreement is only good to about 20%, the 
approximation preserves the general trends pointed to by the detailed 
models. 

Figure 1 presents the results of the approximation for w - - 4 and 
E Q m 0.5 x 10 51 ergs in the form of a contour plot. The vertical axis 
is log T s , the post-shock temperature. The horizontal axis is log n^, 
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where is the normalizing density at R Q " Rj_. It is a constant over 
the time development of a single blast wave, so a given blast wave 
will evolve straight down in this diagram. The contour lines are 
those of the C-band approximation of equation (2.1), and the plus 
signs are the various detailed runs, labeled with their run numbers. 
This entire diagram has E Q fixed at 0.5 x lO"’*' ergs. The kinks in the 
contour lines around log T g « 6.0 show the boundaries between the 
various approximation regimes (equilibrium or non-equilibrium cooling, 
and the cold correction exp(l-l/Tg) for Tg < 1). In reality the 
transitions would be smoother, but it is useful to be able to see on 
the diagrams where the changes in the approximation occur. 

Model 17 falls in a regime which is interesting for several 
reasons, so it will be described in detail here. The ionization state 
of the gas rises rapidly just behind the shock; oxygen ions will be 
given as an example. Moving inward from the edge, the fraction 
n(0 VII)/n(0) rises to essentially unity by R/R„ = 0.96. 0 VII 
remains the dominant species to R/R “ 0.80, including more than 99 % 
of the mass. Nevertheless, there is enough 0 VI in the region 
0.96 < R/Rg <1.0 to produce a column density of 3.6 x 10^ cm This 
0 VJ. will have a temperature of approximately T g (0.56 x 10^ K In this 
case), a velocity of roughly 3/4 of the shock velocity, or about 
150 km/s. Jenkins (1978a, b) used the Copernicus Satellite to observe 
0 VI column densities toward 72 early-type stars at various distances 
from 85 pc to “2 kpc. These measurements tend to show detections with 
column densities N(0 VI) ~ (1-2) x 10^ c,m~^, with velocities near 
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zero (with a star-to-Btar dispersion of 26 km s”*) and thermal widtho 
Indicative of temperatures around 3 x 1G^ K. Some of the stars Bhow 
only upper limits of roughly this magnitude. There is some difficulty 
fitting stellar continua near these lines, so that very broad profiles 
might have been missed (Jenkins, 1984). The column density 
discrepancy will be discussed in a moment; the velocity problem 
remains unresolved. 

Of the power radiated in the B and C bands, roughly 2/3 is 
emitted in a dozen lines of Si VII, Si VIII, Si IX, Mg VIII and Fe X. 
Hence the abundances assumed for these heavy elements are crucial to 
the calculations, as is the assumption that these elements are in the 
gas phase prior to being shocked. If these elements must be sputtered 
from silicate dust grains, they will probably not reach such high 
ionization stages as quickly, which might also make the emission less 
bright. 

One problem that arises In the calculation of the numerical 
models is that they were assumed to be "non-radlative" (that Is, 
radiation losses were neglected in the dynamical calculations). As 
the temperature becomes lower and lower, however, this assumption 
eventually fails (since gas with log T between 5 and 6 radiates with 
great efficiency [Raymond, Cox, and Smith, 1976]). In practice, the 
radiation rates for the various gas parcels were calculated only at 
the end of the remnant evolution (since we assumed for dynamical 
purposes that no radiation has taken place previously). Hence we are, 



strictly speaking, unable to calculate the amount of energy previously 
radiated by a gas parcel. However, a reasonable approximation (and 
probably an overestimate) to this qunnU':y can be obtained by 
multiplying the final radiation rate by the time since the parcel was 
shocked. This quantity can then be compared to kT, so that some idea 
can be formed concerning which parcels have cooled significantly. The 
resulting quantity typically rises quite sharply and then falls off 
steeply behind the shock, so that even though the calculation is very 
approximate, the acfuwl result (namely the identification of cooled 
parcels) is not very sensitive to the details of the approximation. 

Once the cooled parcels have been identified (and these are, 
where they exist at all, invariably near to the shock and few in 
number), they are excluded from the line integrals used in calculating 
the surface brightness. The runs marked in table 2 have been 
corrected in this manner. 

This approximation probably yields excellent results for the 
C-band, and likely bIbo the B-band, but not for softer X-rays which 
are emitted copiously by the cooler gas nearest the shock front (i.e. 
by gas with log T 5.5), where the approximation is weakest. 

This approach also provides an Interesting insight into the 0 VI 
problem. CA found that the column densities of 0 VI predicted by 
their models were quite insensitive to all input parameters (other 
than the assumed pre-shock ionization state of the gas), and further 
that these were in excess of those reported by Jenkins (1978a, b). 
The present program suffers from the same effect. In the cooler 
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models, however, most of the 0 VI is in the outer shells which are 
identified by the above procedure as having cooled. It is expected 
that substantial recombination will have occurred in these regions 
beyond what our (non-radiative) models indicate. If one thus Ignores 
the 0 VJ. which these models predict in gas parcels that in fact have 
cooled, the column densities of 0 VI are brought down to values of 
(1-2) x lO* 3 cm**'*, comparable to the limits of Jenkins (1978a, b). 

III. Anisotropic Media and loobarichrones 

One obvious feature of the soft X-ray background Is an 
anti-correlation between the B and C bands and neutral hydrogen column 
densities as measured at 21 cm. Sanders et. al. (1977) argue that 
this effect is not due to absorption of the X-rays by Intervening 
neutral gas, since both the best-fit value and lack of energy 
dependence of the cross section fail to match those predicted by 
atomic physics. If the gas were clumped into optically thick clouds, 
an absorption picture can be made to work, but Jahoda et. al. (1985) 
find from their 21 cm observations that the gas is not sufficiently 
clumped to produce the observed anti-correlation with soft X-rays. An 
alternative model Is proposed by Sanders et al. (1977) who point out 
that volume occupied by hot X-ray emitting gas is not occupied by 
neutral absorbing gas. Hence at a given galactic latitude, if in a 
given direction one observes a small 21 cm column density, more room 
is left for X-ray emitting gas, which increases the emission measure 
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nnd hence Che surface brightness of soft X-rays. This model aflijwaes a 
uniform emitting region. 

The present work attempts to put more physical verisimilitude 
into this latter displacement picture. The cavity models suggested a 
scheme for examining in an approximate way, using paired spherically 
symmetric models, the effects of a non-isotropic ambient medium, such 
as one containing large clouds. While the blast wave Interacts only 
with the front surface of the cloud, there is some reason to think (as 
in the displacement picture of Sanders et al.) that if the neutral 
material begins further from the Sun, there remains less room for it 
before the line of sight breaks out of the disk of the Calaxy, and 
hence the column density will be smaller. The models actually examine 
the effects on X-ray flux of the pre-shock density and the distance to 
the shock front. 

Suppose the preexisting cavity has the same power-law density 
dependence in all directions, but with a larger constant of 
proportionality in some directions (toward "clouds") than others. 
Consider two directions A and B, with A being toward a denser region 
of the sky. Then let the ambient density laws be 


n oA 


n A R“ w 


n oB 


n B 


r-W 


( 3 . 1 ) 
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ln the two directions, with n A > n B and the same value of w. A 
schematic diagram of this situation is presented in figure 2. 

The interior temperatures (and hence the sound speed) are very 
high, so if two separate models are to be used to represent the same 
blast wave in these two directions, they share the same center and 
thus must have the same interior pressure. Further, they must have 
had the same interior pressure throughout their histories. The two 
models should thus be characterized by the same age. These two 
requirements suggest the name "Isobarichrones" for two such models: 
they have the same age and interior pressure history. 

In Paper II it was shown that time and edge pressure develop 
approximately as power-laws with normalized radius z: 

t * z < 5 “ w >/ 2 (3.2) 


-3 


Also, the central plateau pressure is always nearly one-third of the 
post-shock pressure, so that p c « z“ 2 . Thus if two models are to have 
both the same age and pressure history, we are constrained to pairs 
with the same value of w, as was assumed above. 
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X. 


Ab .the two models evolve, there are two competing effects on 
their X-ray surface brightness. First, the higher density model A 
tcndi> co be brighter for reasons discussed by CA, basically that more 
material enhances the emission measure. Second, since most of the 
X-ray power is in lines of highly ionized species, the degree of 
ionization is critical to the X-ray band surface brightnesses. Model 
A, since it has a higher ambient density,- also has a lower shock 
velocity, and hence a lower post-shock temperature. The excitations 
then tend to favor lower energy transitions than model B, tending to 
make it less bright at higher energies. 

To make this quantitative, let the ratio of the shock radii be 
given by the parameter X (which is less than one in the present 
notation) : 
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XR 
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(3.3) 


Since (Paper XI) E Q =« 3V s p s e(u)/2 where V g is the volume of the 
remnant, e( a>) is a numerical factor that depends on w (values are 
given in Paper II), and p g and w are the same for the two remnants, 
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Thus the explosion energies will be different for isobarichrone model 


pairs. Qualitatively this means that some of the energy originally 
directed toward the cloud (model A) has washed around the cloud (into 
the model B directions). 

Using equations (27) and (28) of Paper IX (with An A fixed) and 
eliminating and t^ in favor of E Q and n^ (where n Q “ n^z w ], a bit 
of algebra produces 

_j 4— (li 
. 7-3 w 
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for the normalizing density n^, the normalized shock radius z (which 
governs the maturity of the remnant with respect to thermal 
conduction), and the present-day ambient density n Q . The ideal gas 
law then produces the post-shock temperature law 
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(In these models, x fl approaches the familiar value of 4 at large z, In 
fact for 2 slightly greater than unity.) 

Using the approximate C-band count rate function from section II 
above, it is possible to examine the behavior of the surface 
brightness by calculating the derivative dC/d\. This can be done in 
each of the three regimes (non-equilibrium cooling, equilibrium 
cooling, and T g < 10® K). This demonstrates that the first two cases 
produce brighter C-band rates in direction A (the "cloud”) > i. e. that 
C correlates with n Q (and hence presumably with the column density Ny 
as viewed from inside the remnant). Only in the third case 

C 

(T s < 10 g K) does one find the desired anti-correlation. In fact the 
"cold correction" factor exp(l-l/Tg) is responsible for this 
anti-correlation behavior. In other words, this approximation scheme 
predicts that when the post-shock temperature of remnant B is less 
than a million degrees, remnant A will be dimmer. This is borne out 
by several sets of detailed model calculations for isobarichrone sets. 

Figure 3 is a contour plot which illustrates this situation for 
E q jj 5 x 10®^ ergs and w =* -4. The axes are Tgg, the post-shock 
temperature of a remnant chosen here to have C » 200 counts per second 
(direction B), and R SA« the minor shock radius (direction A). 
Applicable detailed model runs are plotted as plus signs and marked 
with their run numbers. Isobarichrone families of models fall along 
horizontal lines in this diagram, and runs 13 and 17 are "B" models 
corresponding to the "A" models 14 (for 13) and 18 and 19 (for 17). 
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The contour lines are the approximate C-band rates from equation 

(2.1) . The 200 count s" 1 contour which runs from upper left to lower 
right is the \ m 1 line, where the radii (and densities) in the two 
directions A and B are equal. To the left of this line, A < 1, so 
RgA < Rg B and n Q ^ > n Q g. These contours show a ridge running from 
upper left to lower right, with a saddle at approximately the location 

i 

of run 17. (The 200 count s contours should cross at the saddle 
with \ - 1 following the more vertical one.) ThuB below T gB - 10^ K, 
we have an anticorrelation between C and n Q (i. e.. a smaller R g ^ 

f. 

gives a smaller C-band rate), while the reverse is true above 10° K. 

Clearly, since all this analysis was done with spherically 
symmetric models, this approximation is sufficiently crude that little 

t 

if anything can be said about the edges of the clouds. 

IV. A Summary: Contour Plots 

As a summary of the behavior of the models we have calculated, 
and in an attempt to constrain the model parameters for the local 
bubble, we present several contour plots of parameter space. Equation’ 

(2.1) presents a manner of estimating the C-band intensity as a 
function of three parameters which (using the relations derived in 
Paper II) can be taken to be R s , n Q , and E Q . Given these three 
parameters, all other dynamic quantities of interest can be 
calculated. In particular, we can calculate the postshock temperature 
T s , the central pressure which the. local interstellar medium might be 
expected to experience, and (using equation [2.2]) the B/C ratio. 
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The all-sky average value for the C-band intensity is 
approximately 200 counts s"^. We have therefore set C - 200 for 
purposes of these plots, picked values for R s and n Q , and iteratively 
solved for E Q . The other parameters then follow directly as above. 
The resulting contour plots are presented for several parameters as 
figure 4 (for u) a -4). 

Each plot has the post-Bhock temperature contours plotted for 
reference (labeled with the temperature in units of 10° K) in dashed 
lines. Figure 4a also has contours of E,^ » Eq/ 10^ 1 ergs. The upper 
right part of these diagrams is the ’'cool" portion, where the cold 
correction exp(l-l/Tg) is used. In the lower left, the 
non-equilibrium cooling approximation was used. The shaded portion of 
these figures represents a fold in the surface: all of the 
approximation regimes coexist in this area. In other words, for some 
choices of R fl and n Q , there are three choices of E Q (and hence T fl and 
the other parameters) which will produce C “ 200, depending on which 
approximate cooling function one chooses. The Tg » 1 contour is In 
the shaded fold zone. 

In the extreme upper right part of these diagrams (Tg < 0.4), the 
edge of the blast wave will have cooled, bo both the analytical 
approximation (equation 2.1) and the models on which it is based will 
fail. 

Several authors (see Paresce 1984 and references therein) have 
observed that the Sun seems to be in a cavity of extent “100 pc In 
most directions. These results come from UV absorption line studies 
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and extinction measurements. We thus expect that the radluB of the 
present blast wave will be of this order. Figure 4a demonstrates that 
an explosion energy of 5 x 10 50 ergs or more is needed for such a 
radius. 

In order for the isobarichrone scheme of section III to produce 
an anticorrelation of C-bund intensity with the neutral hydorgen 
column density N u , we need Tg < 1, which is the upper right part of 
these diagrams. 

Figure 4b shows contours of the central pressure, expressed for 
convenience as p/lOOQk, where k is the Boltzmann constant. It then 
has units of cm” 3 K. This plot tends to force our attention toward 
the upper right, in order to reduce the overpressure needed within the 
local bubble. Various estimates of the mean Interstellar pressure 
produce p/k values of order 2000 cm” 3 K (e. g. Spitzer 1978), and 
while we expect a large overpressure within an active blast wave, the 
value of p/k « 30,000 cm” 3 K required by this plot for a post shock 
temperature of 10^ K seems uncomfortably high. However, such high 
pressures have been proposed before (e.g. Williamson et al. 1974, 
Shapiro and Field 1976). Jenkins (1984) points out that while the Sun 
seems to be embedded in a very local partially ionized medium with 
T = 1-2 x 10^ K and n * 0.1 cm” 3 , so that p/k - 2000 cm” 3 K, the 
extent of this region (3-4 pc in some directions [Bruhweiler and Kondo 
1982]) plus its sound speed make the shock crossing time of order 
10 3 yr, which is comparable to the dynamical timescales of the present 
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bubble models. Thus the discrepancy In the pressures may be an effect 
no more transient than the blast wave itself. 

Figure 4c presents a contour plot of the B/C ratio, and it is 
apparent that this ratio is an excellent temperature diagnostic. The 
local component of the soft X— ray background seems to have a B/C ratio 
of between 0.35 and 0.45 (Nousek et al. 1982), which constrains T to 
the range 0.6 to 1.0 x 10^ K. 

We thus find ourselves constrained to a rather narrow stripe of 
parameter space. If one wishes to fit the B/C band without overly 
stressing the pressure constraints, and one wishes to use the 
isobarichrone anticorrelation scheme of section XI, one is forced to 
the conclusion that T 0 is (.8 ±.2) x 10 6 K. We consider it an 
important success of this model that the anticorrelation scheme works 
in a temperature regime which 1$ allowed by the observed B/C ratio. 
Note also that the acceptable values of n Q represented are in a narrow 
range; we find that for C “ 200 counts s’" 1 (with w « -4, 
80 pc < R 0 < 120 pc, and 0.6 < Tg < 1.0), we require ambient densities 
of (1.2 ± .2) x 10” 2 cm ”3. 

For other values of w the approximation formula (2.1) has not 
been checked as rigorously against detailed models. The diagrams in 
figure 5 are similar to those in figure 4, only for the to = 0 or 
uniform ambient density case, again with C fixed at 200 cps. These 
plots, and similar diagrams for w “ -2, suggest ambient density ranges 
of (0.7 ± .2) x 10“ 2 cm” 3 for w «* 0 and (0.9 ± .2) x 10 -2 cm” 3 for 
W “ “2. The temperature constraints are Independent of w, producing 
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T 8 - (0.8 ± .2) x 10® K, though the isobarlchrone scheme of section 
III will not work as envisioned in the uniform medium (w ■* 0) case. 

Several questions were suggested by the investigation of CA, some 
of which are addressed by the present project. Among these quesions 
are the following: 

(1) CA produced a Bet of SNR parameters (n Q , R s , and t) that fit 
the B and C band data for E Q - 5 x 10 50 ergs. What ranges can these 
parameters take on without substantially Impacting on the quality of 
the fit? How would these ranges change if E Q were allowed to vary? The 
range of parameters has been addressed above. 

(2) Is it possible (for example, by setting off the explosion in 
a pre-existing cavity) to make models that will produce a more 
significant fraction of the observed M-band flux? The present work 
suggests that explosions in cavities are no more productive of M-band 
X-rays than those in uniform ambient media. This still presents no 
serious difficulty to the local bubble models, but will be more 
serious in attempting to model the North Polar Spur region, which is 
potentially such a large bubble viewed from the outside. 

(3) Are there parameter choices for which the X-rays could be 
produced without at the same time generating a large local 0 VI 
component? It is still true that the active blast wave models 
inherently generate N(0 VI) ~ 2 to 5 x 10^ cm” This component Is 
hot, has large radial velocity, and is found only along a thin zone 
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close to the shock. If, however, the edge of the blest wave has 
cooled, it seems likely that It will have also recombined to stages 
below 0 VI. This appears promising, though more detailed models which 
incorporate radiation effects in a more natural and satisfactory 
manner would shed more light on the question. 
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Figure Captions 

Fig. 1. The analytical approximation of equation (2.1) aa a contour plot of 
C-band aurfacc brightncoa in counta s” 1 , The vertical axis iB log 
T 8 , the poot-shock temperature. The horizontal axio iB log n^, the 
normalizing density of Paper II, which in a constant over the life 
of q given remnant. E Q - 0.5 x LO 51 ergo, and w ■ -4 
(i.e. n Q « R^) . Also plotted arc the positions of applicable model 
runs, with their run numbers. A remnant evolves straight down on 
this diagram. 

Fig. 2. A schematic representation of the cavity discussed in section III. 
Direction A represents a pre-existing cloud. 

Fig. 3. A contour plot of the C-band intensity in direction A (see Fig. 2). 

Axes are log T oB , the post-shock temperature in direction B (which 
is taken to have C - 200 counts s'* 1 ), and R qA , the shock radius in 
direction A. Applicable runs are plotted. Isobarichrone models 
fall on horizontal lines in this diagram. E Q * 0.5 x lO^ 1 ergs, and 
w “ “4. 

Fig. 4* For C •* 200 counts s - * 1 , contours of the post-shock temperature Tg 
(dashed lines, all plots) in units of 10 6 K are presented. The 
shaded region represents a fold in the surface. Axes are 6hock 
radius and current ambient density. a> « —4. 

a) Explosion energy E 51 - E o /10 51 ergs in solid lines, b) Pressure 
p/lOOOk in solid lines, c) Boron to Carbon band ratio B/C in solid 


lines. 
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Fig. 5. Same ns Fig. A, only for w - 0. C - 200 counts s’* 1 os before. 
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Table 1 

DYNAMICAL PROPERTIES AND N(0 VI) 


Run 

R s n o 

pc 10“3 cm - "* 

T 

A s 

10 6 K 

HanMRaaHUMai 

t 

10 5 yr 

inpinaav 

Z 

iMaaBaaaaa 

x s 

■ ■NBNBaaaBi 

E o 

10 51 erg 

aBMBBBBBBBBH 

N(0 VI) 
lol 3 cm" 2 



Time 

Sequence A (u ■ -4) 




5 

83.9 

2.48 

4.33 

0.36 

0.800 

2.782 

0.500 

• * • 

10 

94.0 

3.90 

1.91 

0.60 

0.896 

2.858 

0.500 

3.1 

1 

104.6 

6.00 

0.77 

0.98 

0.998 

3.353 

0.500 

• « • 

11 

110.0 

7.32 

0.48 

1.24 

1.049 

3.769 

0.500 

4.9 

9 

120.0 

10.37 

0.25 

1.87 

1.145 

3.966 

0.500 

6.7 

6 

136.0 

17.10 

0.10 

3.35 

1.298 

3.997 

0.500 

• • • 



Time 

Sequence B ( co => -4) 




8 

58.0 

5.86 

5.56 

0.22 

0.800 

2.782 

0.500 

• ■ • 

12 

62.1 

7.70 

3.41 

0.30 

0.859 

2.809 

0.500 

2.8 

7 

96.1 

44 . 30 

0.11 

2.27 

1.327 

4.000 

0.500 

• • • 



An Isobarlchrone Pair 

( w " “ 

•4) 



13 

75.0 

11.9 

1.08 

0.602 

0.984 

3.238 

0.500 

2.1 

14 

56.3 

21.2 

0.51 

0.613 

1.052 

3.831 

0.211 

4.2 



An Isobarlchrone Triplet ( w 

a -4) 



17 

83.7 

14.1 

0.56 

0.874 

1.053 

3.789 

0.500 

3.6 

18 

62.8 

25.1 

0.30 

0.888 

1.135 

3.960 

0.211 

2.8 s 

19 

54.4 

33.4 

0.23 

0.894 

1.179 

3.981 

0.137 

1.0 s 



Miscellaneous Runs ( m 

- -4) 




15 

64.5 

6.13 

1.56 

0.450 

0.918 

2.906 

0.211 

3.2 

16 

90.0 

15.96 

3.86 

1.125 

1.102 

3.927 

0.500 

6.7 



A Uniform Ambient Density Case ( w 

- 0) 



20 

55.7 

5.00 

5.11 

0.382 

0.519 

3.255 

0.500 

0.26 


a These figures exclude cooled gas near the shock. 
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Table 2 
BAND RATES 


HMdMMRHMMaM 

Run 

Detailed Models 
Be B C 

M 

laMBBtingiUHMaaonNaiBBiiiBRBHiiHPaitnHtiMnu 

Analytical Approximation 
B C % B % C 



Time Sequence A 





5 

1.1 

5.5 

14.0 

0.90 

7.2 

17.1 

30.3 

22.4 

10 

1.8 

11.4 

31.7 

1,12 

12.1 

34.3 

6.1 

8.1 

1 

9.7 

31.8 

57,2 

1.15 

23.1 

54.7 

-27.3 

-4.4 

11 

6.4 

25.7 

50.1 

0.54 

29.0 

58.8 

12.7 

17.4 

9 

, 8.4 

19.4 

32.1 

0.36 

15.9 

26.3 

-18.1 

-18.1 

6 

3.4 

6.9 

12.1 

0.17 







Time Sequence B 





8 

12.5 

17.5 

55.2 

7.10 

20.1 

69.6 

15.1 

26.1 

12, 

4.3 

33.0 

113.6 

11.80 

28.0 

101.7 

-15.0 

-10.5 

7 

17.4 

32.2 

56.5 

1.38 







An 

Isobarichrone Pair 




13 

7.8 

66.3 

196.9 

5.60 

56.4 

158.9 

-14.9 

-19.3 

14 

30.0 

124.1 

254.2 

3.95 

130.1 

277.4 

4.8 

9.1 



An 

Isobarichrone Triplet 




17 

18.8 

87.7 

189.2 

3.15 

91.9 

204.2 

4.8 

7.9 

18 

29. 2 a 

77. 3 a 

137. 7 a 

2. ll a 

88.2 

149.2 

14.1 

8.4 

19 

13. 9 a 

51. 0 a 

102. 4 a 

1.67 a 

50.9 

85.4 

-0.2 

-16.6 



Miscellaneous Runs 





15 

4.1 

23.6 

59.3 

1.40 

19.4 

53.0 

-18.0 

-10.7 

16 

25.9 

76.1 

138.2 

2.13 

83.5 

157.1 

9.8 

13.7 



A Uniform Ambient Density Case 



. 20 

2.4 

34.8 

136.1 

35.70 

27.0 

108.7 

-22.3 

-20.2 


a These figures exclude cooled gas parcels near the shock. 
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